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ABSTRACT 

We use Wbody/smoothed particle hydrodynamics simulations of encounters between an early-type 
galaxy (ETC) and a late-type galaxy (LTG) to study the effects of hot halo gas on the evolution for 
a case with the mass ratio of the ETC to LTG of 2:1 and the closest approach distance of ^100 kpc. 

We find that the dynamics of the cold disk gas in the tidal bridge and the amount of the newly formed 
stars depend strongly on the existence of a gas halo. In the run of interacting galaxies not having a 
hot gas halo, the gas and stars accreted into the ETG do not include newly formed stars. However, 
in the run using the ETG with a gas halo and the LTG without a gas halo, a shock forms along the 
disk gas tidal bridge and induces star formation near the closest approach. The shock front is parallel 
to a channel along which the cold gas flows toward the center of the ETG. As a result, the ETG can 
accrete star-forming cold gas and newly born stars at and near its center. When both galaxies have 
hot gas halos, a shock is formed between the two gas halos somewhat before the closest approach. 

The shock hinders the growth of the cold gas bridge to the ETG and also ionizes it. Only some of the 
disk stars transfer through the stellar bridge. We conclude that the hot halo gas can give significant 
hydrodynamic effects during distant encounters. 

Subject headings: hydrodynamics — methods: numerical — galaxies: evolution — galaxies: halos — 
galaxies: interactions — galaxies: star formation 


1. INTRODUCTION 


Interacting galaxy pairs have been actively studied 
both obs ervationally and numerically over the past few 
decades (iToomre & Toomre||I972l iBarnes & Hernquist 


1992 1 btruck||20U6| |bmith et al.||2UU7[ among many oth¬ 
ers). While the effects of galaxy mergers have been more 


galaxy merger s 
intensely studied ( e.g., [Spr ingel et a l. z 
Genel et al. 2QQ8j biaao et al. ' 2UU9p , the in 


2005 

Gox 

2008: 

the inl 

beresi 

: in a 


whole variety of interactions, including flybys and dis 


only recently (e.g.. 

Park et al.|2008 Hwang & Park||2009l 

Park & Ghoi 

2009 

Park & Hwang|2009t |bmha & Holley- 

Bockelmann 

[2012 

L’Huillier et al. |20I5p. Eor example. 

Park et al. 

pDOS 

) presented a unified picture on the 


distant interactions and m ergers by analyzi ng the Sloan 
Digital Sky Survey (SDSS;|York et al.||2000|) samples. In 
particular, they found the strong dependence of galaxy 
morphology on the local environment set up by the near¬ 
est neighbor galaxy. It was claimed that hydrodynamic 
effects of cold gas as well as hot diffuse halo gas play an 
important role in the evolution of galaxy color and star 
formation activity. It was also shown that galaxy—galaxy 
encounters are very important in a cluster environment in 
star formation quenching through hydrodynamic interac¬ 
tions of late-type galaxies (LTGs) with cluster early-type 
galaxies (ETGs) |Park & Hwang||2009|). 

So far, however, the overall eftects of hot halo gas 
on galaxy evolution during various interactions have not 
been well understood. They can be tackled with the use 
of hydro dynamical models that includ e both cold and 
hot gas. Recently, Moster et al. (2011) studied the role 
of hot gas in major mergers of disk galaxies, including 


for the first time a gradually cooling hot gas halo in the 
galactic-scale simulations. They showed that the hot gas 
component affects the star formation rates (SERs), the 
efficiency of the starburst, and the kinematics and inter¬ 
nal str ucture of the merg er remnants. In a subsequent 
study ( Moster et al.|2QI2 ), they also examined the effects 
of hot gas on disk thickening in minor mergers. They 
found that a cooling hot halo acts to reduce thickening 
by re-forming a massive thin disk in the post-merger. 
These studies have addressed that the results of inter¬ 
actions can be significantly affected by a hot gas halo 
component. 

In order to further investigate the effects of hot halo 
gas in a wide range of interactions, we have been per¬ 
forming a series of Wbody/hydrodynamic simulations, 
including a cooling diffuse hot gas halo in our galaxy 
models. W e presented the resul ts of the first set of sim¬ 
ulations in Hwang et al. (2013) (hereafter Paper I). In 
that work, we tested different density profiles for the gas 
halo with varying mass and initial spin in our Milky Way 
Galaxy—like models and examined how the adopted gas 
halos create differences in the evolution of the isolated 
galaxies, particularly in the gas dissipation and the star 
formation activity in the disk. 

In the current work, we perform a set of simulations 
of interactions between an ETG and an LTG . This work 
is mot ivated by the observational study by [Park etld 
(2008), where some interesting examples of ETGs hav- 


ing a late-type neighbor were presented, such as a blue 
elliptical-looking galaxy with strong star formation at its 
center, and a red galaxy having spiral-looking arms (see 
Eigure 4 of their paper). In this work we aim to provide 
a set of simulated interacting galaxy systems that can 
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be compared with observations in the follow-up studies. 
Our study focuses on, for the first time, the effects of 
hot halo gas during distant encounters on the changes 
of galactic properties. Specifically, we intend to inves¬ 
tigate whether and how the presence of a hot gas halo 
component can cause differences in the star formation ac¬ 
tivities, the development and location of the tidal bridge, 
and the amount and dynamics of the material transfer¬ 
ring through the bridge to the ETG, as in the examples 
mentioned above, with a fixed set of orbital parameters. 

We will continue to study more cases of various inter¬ 
actions in different environments, emphasizing the role of 
hot halo gas in follow-up studies. For example, we plan 
to work on galaxy interactions in a rich galaxy cluster en¬ 
vironment and investigate star formation quenching and 
the morphology transformation of LTGs influenced by 
the hot halo gas of neighboring ETGs, as well as the hot 
cluster gas. 

This paper is organized as follows. We describe our 
initial galaxy models and the model parameters in Sec¬ 
tion 2, and the simulation code used for this work in Sec¬ 
tion 3. In Section 4, we present our simulation results 
of the distant encounters, explaining the orbital parame¬ 
ters, the evolution, and the star formation activities. Fi¬ 
nally, we summarize and discuss our results in Section 5. 


2. INITIAL GALAXY MODELS 

We use four different galaxy models for this work, 
among which two are LTG models L and LH, and the 
other two are ETG models E and EH (Table 1). We gen- 
erate the mod els by using the ZENO software package 
(Barnes 2011). ZENO allows one to construct multi- 
ple spherical and disk components of a galaxy system in 
mutual equilibrium with user-specified density profiles in 
collisionless or gaseous form. 


2.1. Density Profiles 

We first present the specific form of the density profile 
adopt ed for each component o f our galaxy models (re¬ 
fer to [Barnes fc Hibbard||2QQ^ for more details; see also 
Paper TJ^ 

Both star and gas disks in the LTG models follow an 
exponential surface density profile and a sech^ vertical 
profile: 




where the subscript “c” stands for either “s” for the star 
disk component or “g” for the gas disk component; R 
represents the cylindrical radius (x^ M^c is the 

mass of the star or gas disk; Udc and Zdc are the radial 
scale length and the vertical scale height of the disk, re¬ 
spectively; and Udff is chosen to be 2.5 times larger than 
ads (Table 1). 

The dark matter (DM) halo in all of our models has 
has a Na varro—Frenk—White (NEW) profile (Navarro et 
al. 1996). Because the cumulative mass of the profile 
diverges at large radii, we apply an exponential taper at 


radii larger than 5hd- 


Phd{r) 


Mhdidhd) 1 


47r(ln(2) — 



i) r(r + ahd)^ 
2 

g-2/3(r/6hd-l) 


for r < 6hd , 


for r > 6hd , 


( 2 ) 

where Uhd is the radial scale length of the DM halo, and 
Tfhd(<^hd) is the mass within ahd- 
Two of our models (LH and EH) possess a nonsingular 
isothermal gas halo with a taper: 


Phg{r) 


/norm-^hg 1 g-(’’/f'hg)^ 

^hg 


( 3 ) 


where Uhg is the radius of the inner core with a constant 
density, 6hg is the radius of taper, and Mhg is the mass 
of the gas halo. The isothermal profile is chosen because 
it is similar to the fi profile that is commonly used to fit 
the hot gas observed in galaxy clusters. 

The mass dens ity of the bulge in all models follows a 
Hernquistj (|1990|) profile with truncation at large radii 
from 6b (Table 1): 


Ph{r) 


Qb-^b 1 
27r r(ab + r)^ 



for r < 6b , 


for r > 6b , 


( 4 ) 


where Ub is the radial scale length of the bulge; Mb is 
the bulge mass. The bulge consists of stars only. 

With the above density profiles, the Wbody realiza¬ 
tions of our galaxy models are made by using the ZENO 
programs. We have run each model in isolation to check 
the stability before performing simulations of distant en¬ 
counters. In Appendix A, we present the evolution of 
model LH, which possesses both a gas disk and a gas 
halo. 


2.2. LTG Models 

The total mass of each LTG model is set to Mtot = 
126 X 10^^ Mq (Table 1), similar to that of the Milky 
Way Galaxy. 

Model L consists of the four components: stellar disk, 
gas disk, DM halo, and bulge. The mass of each com¬ 
ponent in units of 10^^ Mq is Mds = 4.4, Mdg = 0.6, 
Mhd = 120, and Mb = 1, respectively (Table 1). The to¬ 
tal disk mass is thus Md = Mds + Mdg = 5.0 x 10^^ Mq. 
The gas fraction in the disk, defined as /dg = Mdg/Md, 
is 0.12. The number of star or gas particles distributed 
in each component is Ms = 204,800, Mg = 81,920, 
Md = 1, 024, 000, and M = 45, 056. 

Model LH has an additional component of a gaseous 
halo. The gas halo component has Mhg = 1.2 x 10^0 Mo 
and Mg = 163, 840. Because the DM halo mass is set 
to Mhd = 118.8 X 10^^ Mq, the total halo mass is Mh = 
120 X 10^^ Mo, equal to that of model L. The halo gas 
fraction is /hg = Mhg/Mh = 0.01. The gas halo is twice 
more massive than the gas disk in this model. The mass 
of a single gas particle (either in the disk or in the halo) 
is the same at the initial time. 

Both star and gas disks in the LTG models have an 
exponential surface density profile and a sech^ vertical 
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Table 1 

Initial galaxy models 


Model components 

Model parameters 

Model L 

Model LH 

Model E 

Model EH 


(IQlO Mq) 

126 

126 

252 

252 


/dg (disk gas fraction)^ 

0.12 

0.12 




/hg (halo gas fraction)^ 


0.01 


0.01 

Star disk: 

Disk model 


Exponential 

Exponential 



Uds (kpc) 

Star disk scale length 

3.5 

3.5 



(kpc) 

Vertical disk scale height 

0.35 

0.35 



^ds (kpc) 

Outer disk cutoff radius 

42 

42 



Md, (IQio Mq) 

Total star disk mass 

4.4 

4.4 



A^ds 

Number of particles 

204 800 

204 800 



mds (10® Mq) 

Mass of individual particles 

2.15 

2.15 



Cds (kpc) 

Gravitational softening length 

0.07 

0.07 



Gas disk: 

Disk model 


Exponential 

Exponential 



ttdg (kpc) 

Gas disk scale length 

8.75 

8.75 



Zdg (kpc) 

Vertical disk scale height 

0.35 

0.35 



6dg (kpc) 

Outer disk cutoff radius 

105 

105 



Mdg (IQio Mo) 

Total gas disk mass 

0.6 

0.6 



iVdg 

Number of particles 

81 920 

81 920 



mdg (10^ Mq) 

Mass of individual particles 

0.73 

0.73 



Cdg (kpc) 

Gravitational softening length 

0.04 

0.04 



DM halo: 

Halo model 


NEW 

NEW 

NEW 

NEW 

tthd (kpc) 

DM halo scale length 

21 

21 

21 

21 

^hd (kpc) 

Tapering radius 

84 

84 

84 

84 

Mhd(ahd) (10^° Mo) 

Mass within radius ahd 

12.35 

12.23 

24.71 

24.46 

Mhd(oo) = Mhd (IQio Mo) 

Total DM halo mass 

120 

118.8 

240 

237.6 

iVhd 

Number of particles 

1 024 000 

1 024 000 

2 048 000 

2 048 000 

mhd (10^ Mo) 

Mass of individual particles 

11.7 

11.6 

11.7 

11.6 

Chd (kpc) 

Gravitational softening length 

0.16 

0.16 

0.16 

0.16 

Gas halo: 

Halo model 



Isothermal 


Isothermal 

dhg (kpc) 

Gore radius 


10.5 


10.5 

6hg (kpc) 

Tapering radius 


420 


420 

Mhg (IQio Mo) 

Total gas halo mass 


1.2 


2.4 

^hg 

Number of particles 


163 840 


327 680 

mhg (10^ Mo) 

Mass of individual particles 


0.73 


0.73 

Chg (kpc) 

Gravitational softening length 


0.04 


0.04 


Bulge: 


Bulge model 


Hernquist 

Hernquist 

Hernquist 

Hernquist 

ttb (kpc) 

Bulge scale length 

0.7 

0.7 

2.8 

2.8 

bh (kpc) 

Truncation radius 

140 

140 

280 

280 

Mb (lOio Mq) 

Total bulge mass 

1 

1 

12 

12 

Nb 

Number of particles 

45 056 

45 056 

540 672 

540 672 

mb (1O®M0) 

Mass of individual particles 

2.22 

2.22 

2.22 

2.22 

Cb (kpc) 

Gravitational softening length 

0.07 

0.07 

0.07 

0.07 


^ Mtot — — ^ds + ^dg + ^hd + ^hg + 

^ /dg = Mdg/Md 
c = Mhg/Mh 


profile (Equation (1)). The radial and vertical scale 
lengths of the star disk are set to a^s = 3.5 kpc and 
^ds = 0-35 kpc, respectively. Those of the gas disk 
are set to ang = 2. 5 x a^s = 8.75 kpc and z^g = 
(Moster et al. 2011). The mass distribution of the DM 
halo follows the hJEW profile (Equation (2)). Its radial 
scale is Ubd = 21 kpc, which is adopted from |McMillan 


scale IS ahd = zd Kpc, wnicn is aaoptea iroiii |ivicivinian 
fc Dehn^ (2007). Eor the bulge, the Hernquist pro- 


file (Equation (d)) is chosen with the length scale of 
Ub = 0.7 kpc (jMcMillan & Dehnen||2007|). The gas halo 
included in model LH follows the isothermal density pro¬ 
file (Equation (3)) with the core radius of Uhg = 10.5 kpc. 


The disk gas particles in models L and LH are set to 
rotate in the clockwise direction with the local circular 
velocities without any motion in the vertical direction at 
the initial time. Eor the disk star particles in both mod¬ 
els, velocity dispersions are added to the lo c al cir cular 
velocities, as described in Barnes & Hibbard (2009) ^ 
The temperatures of the disk gas particles in botn mod¬ 
els are set to T = 10, 000 K. Eor the halo gas particles in 


^ In Paper I, the dispersions were added to the disk star particles 
in the models of type DH, but were not applied in the models of 
type D. 












































4 


Hwang & Park 


Table 2 

Three runs of distant encounters 




ETG 

LTG 



1st closest approach 

ETG-LTG 

xo,yo,zo 

(kpc) 

'^xO: '^yO^ '^zO 

(km s-l) 

xo,yo,zo 

(kpc) 

XxOi XyQ., VzO 

(km s-l) 

dr 

(kpc) 

dv 

(km s-i) 

t 

(Gyr) 

E-L 

0, 0, 0 

0, 0, 0 

-1200, 265, 0 

200, 0, 0 

94.3 

495.6 

4.75 

EH-L 

0, 0, 0 

0, 0, 0 

-1200, 265, 0 

200, 0, 0 

94.4 

494.8 

4.75 

EH-LH 

0, 0, 0 

0, 0, 0 

-1200, 265, 0 

200, 0, 0 

94.8 

495.5 

4.75 


model LH, the initial temperatures are determined in ac¬ 
cordance with the hydrodynamic equilibrium of the halo 
gas0 

2.3. ETG Models 

The ETG models E and EH do not have a disk com¬ 
ponent. The total mass of each ETG model is set to 
252 X 10^^ Mq (Table 1), twice more massive than that 
of the LTG models. 

Model E consists of two components, a DM halo 
with Mhd = 240 X 10^^ Mq and a bulge with Mb = 
12 X 10^^ Mq (Table 1). The number of particles dis¬ 
tributed in each component is A^hd = 2,048,000 and 
N\y = 540, 672. The masses of a particle in the DM halo 
and in the bulge of model E are equal to those of model L. 

Model EH has a gas halo in addition to a DM halo and 
a bulge. The gas halo mass is Mhg = 2.4 x 10^^ Mq with 
the halo gas fraction /hg = 0.01. The number of halo gas 
particles is A^hg = 327,680. The mass of the DM halo 
is Mhd = 237.6 X 10^^ Mq. The total (DM + gas) halo 
mass and the bulge mass are equal to those of model E. 

In the ETG models, the DM halo and the bulge follow 
the NEW profile (Equation (2)) and the Hernquist profile 
(Equation (4)), respectively, as in the LTG models. The 
scale length of the DM halo (uhd) and the core radius of 
the gas halo (ahg) are equal to those of model LH (Table 
1). The scale length of the bulge component (ub) is set 
to 2.8 kpc. 


3. SIMULATION CODE AND THE PARAMETERS 

To perform simulations of interactions between the 
galaxy models, we use an early version of the GADGET - 
3 TreeSPH code (originally described in Springel| 2005). 
Gravitati onal force is calculate d by means of the tree al- 
gorithm (Barnes & 


1986), and the hydrodynamic 


force is computed with the smoothed particle hydrody¬ 
namics (^PH)jn£ttmdJiithnentropy conservative formu¬ 
lation (Springel & Hernquist 2002). The code includes 
radiative cooling and heating tor r primord ial mixture 
of hy drogen and helium by photoionization (Katz et al 
1996). Star formation and supernova feedback in the 
interstellar medium (ISM) are implemented using the ef - 
fective multiphase model of|Springel & Hernquist|(|2003|). 
A thermal instability operates for gas at high density 
over a threshold value of pth, and the ISM is treated 
as a statistical mixture of cold clouds and ambient hot 


medium. Star formation occurs in dense regions consum¬ 
ing the cold clouds. The co nsumption times cale is chosen 
to match the observations (|Kennicutt|1998|) . Among the 
newly formed stars, the rnass fraction of massive stars 
(> 8 Mq) is given by the parameter /3. The massive 
stars die instantly as supernovae releasing energy as heat 
to the ambient diffuse medium. Some cold clouds are also 
evaporated inside the supernova bubbles returning ma¬ 
terial to the ambient gas. We do not include the galactic 
winds associated with star formation in this work. We 
also do not consider active galactic nucleus feedbacks. 

Eor the parameters related to star formation and the 
feedbacks, we use standard values of the multiphase 
model in all of our simulations: the star formation 
timescale is tg = 2.1 Gyr, and the mass fraction of mas¬ 
sive stars is p = 0.1. The “supernova temperature” is 
^SN = 10^ K, and the temperature of cold clouds is 
Tc = 1000 K. The parameter value for supernova evapo¬ 
ration is Aq = 1000. In the code, the mass of a gas par¬ 
ticle involving star formation is chosen to be divided in 
two: a new star particle has half of the original gas mass, 
and the mass of the gas particle decreases by half. Thus, 
the total number of particles in each model increases as 
star formation occurs. 

We set the gravitational softening length for the DM 
halo particles to ehd = 0.16 kpc in all simulations. The 
value is about an order of magnitude smaller than the av¬ 
erage distance between the DM particles in our models. 
Then we set the maximum acceleration experienced by 
a particle (/ oc m/e^, with the initial mass) to be equal 
in each component. Eor example, in the simulations of 
model LH in isolation and in interaction with model EH, 
the softening lengths for the gas (both halo gas and disk 
gas), disk star, and bulge particles are determined ac¬ 
cordingly to be 0.04, 0.07, and 0.07 kpc, respectively. 
The softening lengths for disk gas and disk star particles 
are comparable to the average distances between parti¬ 
cles in the components. 

To achieve better accuracy in the force computation, 
we chose the value of the parameter a, which controls 


the a ccuracy of the relative cell-opening criteria (Springel 


2005), to be 0.0005 in all simulations, which is about 


an order of magnitude smaller than the values popularly 
used in galactic-scale simulations. We set the opening 
angle to 6> = 0.3, which is used for the first force compu¬ 
tation. 


^ The temperature distribution of the halo gas particles of 
model LH is statistically the same as that of model DHi in Paper 1. 
Compared with the models in Paper I, we use a larger number of 
particles for this work and accordingly adopt smaller gravitational 
softening lengths in all components. Except for the particle num¬ 
bers and the dispersions of disk stars, all of the other key model 
parameters, such as scale length, truncation/tapering radius, and 
core radius, remain unchanged. 


4. RESULTS OF DISTANT ENCOUNTERS 
4.1. Orbital Parameters and Trajeetories 

We run three simulations of distant encounters between 
the ETG and the LTG to examine the effects of the hot 
gas halo. They are the E-L, EH-L, and EH-LH runs. 
Their initial conditions are listed in Table 2. In all runs. 
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Figure 1. Left: distance (dr^ solid line) and velocity (dv, dashed line) of the late-type galaxy relative to the early-type galaxy with respect 
to time in run E-L. Here, dr is measured between the center of mass positions of the bulges of the two galaxies, and dv is the average 
velocity of the bulge particles of the late-type galaxy relative to that of the early-type galaxy. The closest approach between the two 
galaxies occurs a.t t = 4.75 Gyr and is marked with the vertical dotted line. Right: orbital trajectories of the two galaxies in run E-L seen 
in the x—y plane. The late-type galaxy L starts from the left going to the right (thick solid line) and the early-type galaxy E starts from 
the origin and moves upward (thin solid line), as indicated with the arrowheads. The positions of the galaxies at t = 4, 5, 6, and 7 Gyr are 
marked with circles, crosses, triangles, and squares on the trajectories, respectively. The closest approach between the galaxies is marked 
with the dotted line. 


the mass ratio of the ETG to LTG is 2:1 (Table 1). (We 
have also considered the E-LH case. This run is presented 
in Appendix B.) 

In each run, the ETG is initially placed at the 
origin and the LTG is positioned at {xo^yo^zo) = 
( — 1200 kpc, 265 kpc, 0) with the initial relative velocity 
of {vxo,Vyo,Vzo) = (200 km s“^,0,0) toward the ETG. 
The initial y position (initial impact parameter) of the 
LTG is chosen so that the closest approach distance be¬ 
comes about 100 kpc. (Throughout this paper, we con¬ 
sider the position of the center of mass of the bulge of 
each galaxy as the center position of the total galactic 
system.) The disks of all LTGs are initially set in the 
x—y plane without tilt and have clockwise directional 
spin when viewed from the positive z axis. 

The three simulations are run for over 8 Gyr. In the 
left panel of Eigure 1, the distance between the LTG and 
the ETG, dr, and the relative velocity of the LTG, dr, 
in run E-L are plotted as a function of time. The two 
galaxies encounter most closely at t = 4.75 Gyr with 
dr = 94 kpc. (In this paper, time t represents the time 
elapsed since the start of the simulation, if not otherwise 
specified. The projected distance between the galaxies 
onto the x—y plane is close to the three-dimensional dis¬ 
tance dr throughout the run because the galaxies move 
mostly in the x—y plane.) In the right panel of Eigure I, 
the orbital trajectories of the late-type (thick solid line) 
and the early-type (thin solid line) galaxies in run E-L 
are presented. The positions of the two galaxies at t = 
4, 5, 6, and 7 Gyr are marked on their trajectories with 
circles, crosses, triangles, and squares, respectively. 

As for runs EH-L and EH-LH, the orbital trajectories 
of the galaxies are almost the same as those in run E-L 
(see Table 2 for the closest approach events). 

4.2. Evolution of Run E-L 

Eigure 2 (top two rows) shows the snapshots taken at 
t = 4, 5, 6, and 7 Gyr. The closest approach is at t = 
4.75 Gyr. The time measured from the closest approach 
(written in the parentheses) and the three-dimensional 


separation between the two galaxies, dr, at each time 
are given in the lower panels. 

The second snapshot shows early development of a 
tidal bridge and a tail out of the stellar and the gas disks 
of L shortly after the closest approach. The third snap¬ 
shot shows that non-star-forming disk gas (blue points) 
and “old” disk stars (orange points; the disk stars ex¬ 
isted before the closest approach) are transferred from L 
to E. The cold non-star-forming gas flows through the 
tidal bridge and dynamically follows the motion of the 
accreted old stars. However, none of the green points, 
which represent the stars formed out of the gas since the 
closest approach (hereafter, “young” (disk) stars), nor 
any of the cyan points, which represent the gas particles 
having positive values of SER at that instant, i.e. star¬ 
forming gas, are transferred to E in this run. In the last 
snapshot, the accreted cold gas and the old stars appear 
to orbit around the center of E. (Throughout this paper, 
we refer to the particles by their original setup, if not 
otherwise specified. Eor example, “disk star particles of 
the LTG” mean the particles initially set as the disk star 
particles of the LTG. Similarly, “disk gas particles” mean 
the particles initially set as the disk gas and those have 
not turned into stars by the time.) The total mass of 
the old and young disk stars is shown in Eigure 3 (top 
panel). 

To show the mass transfer from L to E more clearly, we 
display in Eigure 4 the disk particles of L together with 
the velocity vectors in a reference frame centered on the 
center of the ETG projected onto the x—y plane. The 
disk star and gas particles move toward E through the 
bridges in a dynamically similar fashion. To estimate 
the amount of the disk particles transferring to E, we 
use a sphere centered on E with the radius to the point 
where the net local gravitational acceleration becomes 
zero along the line connecting the centers of E and L, as¬ 
suming each galaxy as a point mass. At t = 6 Gyr, some 
of the old disk stars (orange points) and a relatively large 
number of cold non-star-forming disk gas (blue points) 
are transferring to E, seen within the sphere; however. 
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Figure 2. Top two rows: four snapshots of the distribution of the particles in run E-L seen in the x—y plane at t = 4, 5, 6, and 7 Gyr 
(from the first to the fourth columns, respectively). The collisionless (DM and star) particles and the gas particles are displayed separately 
in the upper and the lower panels, respectively. In each panel, the positions of the center of mass of the bulge of the early- and the late-type 
galaxies are marked with a circle and a diamond, respectively. In the lower panels, the three-dimensional distance between the two bulges 
at each time is given; the time measured from the closest approach is also written in the parentheses. The different colors are used to 
distinguish the origin of the particles. In the upper panels, the violet and yellow points represent the particles initially set as the halo DM 
and the bulge stars of the early-type galaxy, respectively. The light-blue points denote the particles set as the halo DM of the late-type 
galaxy. The orange points represent the “old” stars initially set as the disk stars and those added onto the disk out of the gas originally set 
as the disk gas before the closest approach, i.e., t <4.75 Gyr. The green points represent the “young” stars formed out of the initial disk 
gas since the closest approach, t > 4.75 Gyr. (The bulge particles of the late-type galaxy are not displayed for simplicity.) In the lower 
panels, the blue points denote the non-star-forming gas particles, which were set as the disk gas of the late-type galaxy (and remain in gas). 
The star-forming gas particles having positive values of star formation rate at that instant are shown with the cyan points. Middle two 
rows: the same snapshots as the upper two rows, but from run EH-L. The red points in the lower panels represent the particles originally 
set as the halo gas of the early-type galaxy EH. Bottom two rows: the same snapshots as above, but from run EH-LH. The pink points 
in the lower panels represent the particles originally set as the halo gas of the late-type galaxy LH. The green points in the upper panels 
include (at the central part of the disk) a small amount of the young stars formed out of the gas originally set as the halo gas of LH. The 
cyan points in the lower panels also include a few star-forming gas particles that were initially set as the halo gas of LH. 
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Figure 3. The growth of the total mass of the stars on the disk of the LTG together with those transferred to the ETG (if there are any) 
in the three runs. In each panel, the orange shade represents the total mass of the old disk stars, i.e., initial disk stars plus stars formed 
before the closest approach. The diagonal lines are drawn to show the total mass of the initial disk stars, which is 4.4 x 10^® Mq in all 
runs. The closest approach time is marked with the vertical dotted line at t = 4.75 Gyr. The total mass of the young stars that have been 
formed since the closest approach are shown with the green shade. 


none of the young stars (green points) nor any of the 
star-forming disk gas (cyan points) are found within the 
sphere. By t = 7 Gyr, old disk stars and cold disk gas 
have continued to transfer to E. Those materials cap¬ 
tured by E move around the center of E. The star and 
gas bridges begin to fade away. Eigure 5 (top panels) 
presents the total mass of the disk materials enclosed 
within the spheres with respect to time after the closest 
approach (for t > 5.5 Gyr). At t = 6 Gyr, about 0.45 x 
10^ Mq of old disk stars and 3.3 x 10^ Mq of non-star¬ 
forming disk gas are accreted onto E. The amount of the 
transferred cold gas is about 5% of the original cold disk 
gas mass. The accreted materials do not include young 
disk stars or star-forming gas in this run. 

4.3. Evolution of Run EH-L 

As shown in the snapshots in Eigure 2 (third and fourth 
rows), the distribution of the old disk stars (orange points 
in the upper panels) is generally the same as that in the 
E-L run throughout the simulation. However, the young 
disk stars (green points) and the disk gas (blue points) 


show very different configurations. In the first snapshot 
at t = 4 Gyr, the gas disk forms a bow-like front as it 
moves against the halo gas included in EH (red points in 
the lower panels), and some disk gas at the other side is 
stripped. In the next snapshot, at 0.25 Gyr after the clos¬ 
est approach, the gas bridge develops nearly straightly, 
while the stellar bridge consisting of the old disk stars 
bends below the gas bridge. There is a burst of star for¬ 
mation (green points) in the shock-compressed gas along 
the gas bridge right after the closest approach. The offset 
between those young stars and the old stars reflects the 
strength of the shock. The gas bridge, including some 
of star-forming gas (cyan points), continues to extend 
almost radially to the center of EH. The last snapshot 
shows that the accreted old and young stars orbit around 
EH with different orbits. Some young stars and cold gas 
accumulate at the center of EH. It should be noted that, 
in the previous run E-L, no young stars manage to trans¬ 
fer to the ETG (see the top two rows in Eigure 2). The 
middle panel of Eigure 3 shows overall more production 
of young stars (including those transferred to EH) in the 
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Figure 4. Disk materials in run E-L with the velocity vectors (relative to the early-type galaxy) projected onto the x—y plane at t = 
5, 6, and 7 Gyr (from the top to the bottom rows, respectively). The disk stars and gas particles are shown in the left and the right 
columns, respectively. The separations between the two galaxies and the times are shown in the left panels. The times written in the 
parentheses are those measured from the closest approach. In each panel, the colors are used in the same way as in Figure 2 (orange for 
the old disk stars, green for the young disk stars, blue for the non-star-forming disk gas, and cyan for the star-forming gas). The arrows 
represent the transverse velocities of the disk particles in a reference frame centered on the center of the ETG. The lengths of the arrows 
are proportional to the magnitude of the transverse velocities. The center-of-mass positions of the early-type and the late-type galaxies 
are marked with a small circle and a small diamond, respectively. A line connecting between the centers of the two galaxies is drawn. On 
the line, at the position where the net local gravitational acceleration becomes zero, assuming each galaxy as a point mass, another line 
is drawn perpendicular to the line connecting the centers. At later times (after t ~ 5.5 Gyr), the particles lying within the large spheres 
inscribed in the perpendicular line are considered to be captured by the early-type galaxy. 
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Figure 5. Disk materials transferred to the early-type galaxy within the spheres shown in Figures 4, 6, and 8 in runs E-L, EH-L, and 
EH-LH (from the top to the bottom rows, respectively). The left panels show the amount of the disk star particles in units of 10^*^ Mq. 
The amounts of the old disk stars and the young stars formed since the closest approach are displayed separately, with the orange and 
the green solid lines, respectively. The total amount of the old and the young stars is shown with the black dashed line. The right panels 
present the amounts of the non-star-forming disk gas (blue solid line), star-forming disk gas (cyan solid line), and the sum of the two (black 
dashed line) in log scale. In both runs E-L and EH-LH, no young stars or star-forming gas transfer to the early-type neighbor. 


EH-L run than in the E-L run. 

The disk particles are displayed separately in Eigure 6. 
The total masses of the disk gas and disk star particles 
transferred to EH are shown in Eigure 5 (middle row). 
Compared to run E-L, a similar amount of the old disk 
stars but less than half the amount of non-star-forming 
cold gas is transferred to EH within the large spheres 
drawn in Eigure 6. At t = 5 Gyr (top row in Eigure 6), 
the young stars (green points) and the star-forming gas 
(cyan points) are seen over a wider area along the bridges 
compared to those in run E-L; some of them are located 
beyond the line of the zero gravity toward EH. At later 
times, EH is able to accrete the star-forming gas and 
some of the young stars to its center. In addition, some 
non-star-forming gas particles accumulated at the center 
of EH later turn into stars in this run. 

In Eigure 7, we present the local density pgas and the 
temperature T of the gas particles located near the or¬ 
bital plane (within | 2 ;| < 15 kpc) to show how cold disk 
gas is affected by hot halo gas through the hydrodynamic 
interaction. As the gas disk approaches the gas halo (top 
row), the gas density at the leading side of the disk rises, 
and some halo gas piles up along the bow-like front of the 
disk. The piled gas has higher temperatures due to the 
shock heating. Shortly after the closest approach (middle 
row), a channel for the cold gas toward the center of E is 
created. This channel is seen between the inner trailing 
shock wound around EH taken over by the supersonic 
motion of L and the preceding shock that appears due to 


the cold gas tidal bridge. The cold gas can flow though 
this channel and is not blocked by shocks because the 
shocks exist parallel to the cold gas bridge. (We have 
checked that the positions of the shocks in this run and 
in run EH-LH, which strongly determine the subsequent 
evolution, are unchanged in different-resolution runs. We 
present the high-resolution simulations in Appendix G.) 

4.4. Evolution of Run EH-LH 

Eigure 2 (bottom two rows) shows the four snapshots 
from run EH-LH, where both ETG and LTG possess hot 
gas in their surrounding halos (red and pink points, re¬ 
spectively). The old disk stars evolve in the same way 
as in the two previous runs at all times, but the cold 
gas evolves in a unique way, being affected by the gas 
halos. At t = 4 Gyr, a large-scale shock is formed be¬ 
tween the two hot gas halos (bottom-left panel). By this 
time, some disk gas at the outskirt has been ionized and 
scattered into the halo because of the hot halo gas of 
its own galaxy. In addition, the gas disk in this run has 
dissipated more quickly than the other two runs and pro¬ 
duced more stars because of the compression exerted by 
its halo (see Paper I for more details), before the gas 
halo of LH is affected by the collision with the gas halo 
of EH. At later times (t = 5 through 7 Gyr), the gas disk 
does not develop a strong bridge as in the other runs, 
and almost no cold gas is transferred to EHjf] Some halo 

^ The strength of a gas bridge will vary in different settings. 
For example, in the case of closer encounters, a gas bridge might 
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Figure 6. The same as Figure 4, but for run EH-L. 


gas particles of LH cool down to the disk and then turn 
into stars. The total mass of these stars formed out of 
the original halo gas of LH is less than ^ 1 % of that of 
the stars generated out of the original disk gas; they are 
found at the central part of the disk. Figure 3 (bottom 

develop more strongly to the neighboring galaxy. 


panel) shows that the total mass of the young stars in¬ 
creases rapidly near the closest approach but slows down 
after a couple of gigayears. These young stars are all 
found on the disk of LH, not within EH. 

Figure 8 displays the distribution of the disk particles. 
At t = 5 Gyr, the young stars formed since the clos- 
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Figure 7. The local density and the temperature of the gas particles in run EH-L at t = 4, 5, and 6 Gyr. In each panel, the gas particles 
of both early- and late-type galaxies near the orbital plane, with | 2 ;| < 15 kpc, are displayed. The colors of the gas particles represent the 
values of log pg^s in 10^° Mq kpc“^ (l eft column) and l opT in K (right column). The ranges of the color bars in each column are fixed. 
This figure is made by using SPLASH ( [Price Eh J.|2007| >. 


cst approach (green points) appear along the bridge, but 
are less extended compared to run EH-L; the gas bridge 
does not grow beyond the line of zero gravitational ac¬ 
celeration. At t = 6 and 7 Gyr, some old stars (orange 
points) are seen transferred to EH. However, no young 
stars and only a little non-star-forming gas (blue points) 
are captured by EH. The total masses of the star and 
gas particles transferred to EH are presented in Eigure 5 
(bottom row). The amount of the old stars transferred 


to the early-type neighbor is similar to (a little bit more 
than) those in the previous runs. But the total mass of 
the cold gas moved to EH is much less than that in the 
other two runs. 

The local gas density and the gas temperature are 
shown in Eigure 9. The discontinuities of the gas quan¬ 
tities are seen at the shock. As LH turns around EH 
(middle row), the inner trailing shock propagates toward 
the center of EH before the cold gas tidal bridge starts 
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Figure 8. The same as Figures 4 and 6, but for run EH-LH. 

halo gas (bottom row). 

4.5. Star Formation Activities 

We present the SFRs in the three runs in Figure 10. In 
run E-L (top panel), the SFR (blue solid line), which is 
the sum of the SFRs from all gas particles originally set 


to extend toward EH, blocking its motion to the center 
of EH. While the cold disk has been already dissipated 
a lot by this time through star formation, the block¬ 
ing by the shock further hinders the growth of the cold 
gas bridge and ionizes it. The shock is later completely 
wound around EH and makes a spiral pattern in the hot 
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Figure 9. The same as Figure 7, but for run EH-LH. 


as the disk gas, rapidly increases right after the closest 
approach at t = 4.75 Gyr (marked by a vertical dotted 
line). This compares with a gradual decline of SER of 
the corresponding isolated L (black dashed line). The 
value at the first peak of the blue line at t = 4.95 Gyr is 
0.46 Mq yr“^, which is about six times the value of the 
isolated galaxy at the same instant. The excess of the 
SER lasts until t ~ 8 Gyr. 

In run EH-L (middle panel), the SER (blue line) also 
abruptly exceeds that of the isolated galaxy L (black line) 
right after the closest approach. The excess is a little 
greater and lasts longer than that in run E-L. The SER 


at the first peak is 0.5 Mq yr“^ at t = 4.91 Gyr, which 
is about seven times that of the isolated galaxy L. (We 
neglect the SERs of the halo gas particles of the ETG 
model EH because they are very small, although not al¬ 
ways zero, with a maximum value of 0.018 Mq yr“^). 
This implies that the enhanced star formation in this 
run is due to the compression of the gas disk through the 
collision with the gas halo of EH. Although some of the 
disk gas particles are heated and stripped, the compres¬ 
sion by the shock initiated during the collision is mainly 
responsible for the cold gas having a higher SER. Among 
the star-forming gas particles of L, some are captured by 
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Figure 10. Evolution of star formation rate in runs E-L, EH-L, and EH-LH (from the top to the bottom panels, respectively). In each 
panel, the solid blue line represents the sum of the SFRs of all gas particles that initially belong to the late-type galaxy; the dashed black 
line is the sum of the SFRs of all gas particles of the late-type galaxy in isolation. The time of the closest approach is marked with the 
vertical dotted line (red). In the middle panel, the shaded region with cyan color denotes the contribution to the total SFR (shown with the 
blue solid line) from the gas particles transferred to the early-type galaxy (within the spheres shown in Figure 6) since t = 5.5 Gyr. There 
is no cyan shaded region in the top and the bottom panels because none of the star-forming gas particles are captured by the early-type 
neighbor. In the bottom panel, the SFRs of the gas particles originally set as the halo gas of LH are included in both the blue and black 
lines, but the contribution is very small (see the text). 


EH after the encounter. The total SFR of the gas parti¬ 
cles transferred to EH is shown with the cyan shade for 
t > 5.5. Some non-star-forming disk gas transferred to 
EH also turns into star-forming gas particles later, when 
they gather densely enough to satisfy the conditions for 
star formation. These particles lead to the increase of 
the shaded region around t = 1 Gyr. 

In run EH-LH (bottom panel), the SFR (blue line) 
also exceeds that of the isolated galaxy LH (black line) 
but much more greatly than in the other two runs, 
with the SFR of 0.81 Mq yr“^ at the first peak at 
t = 4.87 Gyr. The maximum value of the SFR reaches 
to 1.033 Mq yr“^ at t = 5.2 Gyr. Unlike in the others, 
the excess (over the black line) begins earlier at about t 
= 4 Gyr, somewhat before the closest approach, from 
the onset of the collision between the two gas halos. 
The shock-compressed gas halo of the LTG consequently 


delivers the effects onto the cold disk, raising the SFR 
since t ^ 4 Gyr. The great excess on the SFR near the 
closest approach indicates the strong compression from 
the shock generated by the collision. Gompared with 
the SFR of the isolated L, the SFR of the isolated LH is 
generally higher, mainly due to the compression exerted 
by the gas halo of LH; the direct contribution from the 
original halo gas particles of LH to the total SFR is 
small (at most about 1% or less in the total SFR of LH). 

5. SUMMARY AND DISCUSSION 

Using A^-body/SPH simulations, we studied distant 
encounters between an ETG and an LTG for a case of 
the closest approach distance of about 100 kpc with the 
mass ratio of the ETG to LTG of 2:1. We performed 
three comparison simulations that adopt galaxy models 
having or not having a hot gas halo, keeping all of the 
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Figure 11. Top row: gas particles on the disk at t = 5 Gyr from runs EH-L (left) and EH-LH (right). The particles originally set as the 
disk gas but stripped off or scattered into the halo are not shown. The blue and cyan points represent non-star-forming and star-forming 
gas particles, respectively. The pink points seen in the right panel are the original halo gas particles of LH that cooled down onto the disk 
having a zero star formation rate. The positions of the center of mass of the early- and the late-type galaxies are marked with a circle and 
a diamond, respectively. The horizontal axis is set along the line connecting the centers of the two galaxies. The zero point on the x axis is 
chosen to be the position where the net local gravitational acceleration due to the two galaxies becomes zero. Bottom row: average profiles 
of the temperature T (red dotted line), local density pgas (blue dashed line), and pressure P (magenta solid line) for the gas particles 
around the bridge area, located between the two vertical dashed lines drawn in the top panels. To obtain the average quantities, the gas 
particles shown in the top panels are projected onto the x axis, and then the average is calculated in every 0.5 kpc along the x axis, taking 
N nearest particles into account to the center of each bin. The number of particles N is chosen to be five times the average number of gas 
particles in the bridge within a 0.5 kpc length on the x axis. Each of the three gas quantities in the left and the right panels is displayed 
on the same scale: the temperature scale is shown to the left of the left panel and the local density scale is shown to the right of the right 
panel. The pressure scale is chosen to be fit in the same panel. The gas bridge in run EH-L (left) is more extended than that in run EH-LH 
(right) beyond the vertical solid line at a: = 0. 


other parameters fixed. We summarize our key findings 
in each run as follows, focusing on the mass transfer of 
the disk materials to the ETG and the star formation 
enhancement that is due to the encounter. 

1. E-L Case 

A large amount of cold gas is transferred from L to E, 
but it does not contain star-forming gas. The cold gas 
flows through the gas bridge and dynamically follows the 
motion of accreted old disk stars. Because there is no 
compression of the cold gas by shock, no new stars form 
along the tidal bridge. Thus, the ETG accretes cold gas 
and old disk stars but no young stars nor star-forming gas 
out of the LTG. The accreted cold gas and old disk stars 
orbit around the center of E. The SER of L in run E-L 
increases abruptly right after the closest approach. 

2. EH-L Case 

Near the closest approach, a shock forms along the cold 
disk gas tidal bridge, and new stars form along the 
bridge. The strength of this shock is seen from the off¬ 
set between the tidal bridge that consists of the old disk 
stars and that consisting of the young stars and the cold 
gas. A very important feature of this type of interaction 
is that a channel for the cold gas toward the center of 


EH is created right after the closet approach. This chan¬ 
nel is developed between the inner trailing shock wound 
around EH taken over by the supersonic motion of L and 
the preceding shock that appears due to the cold gas tidal 
bridge. The cold gas can flow almost radially through 
this channel and is not blocked by shocks because the 
shocks exist parallel to the cold gas tidal bridge. The 
cold gas flowing to EH contains some star-forming gas. 
Some of the non-star-forming gas transferred to EH later 
turns into stars as well. As a result, EH can accrete the 
cold disk gas, including some star-forming gas, at its cen¬ 
ter in this particular run, and also accrete the old and 
young disk stars around the center. Some of the young 
stars are found at the center, as well. The SER of L 
in this run also increases rapidly right after the closest 
approach. The increase is greater than that in run E- 
L. Then, the SER decreases most slowly among all three 
runs mainly because the accreted cold gas turns into stars 
later at the center of EH. 

3. EH-LH Case 

As LH approaches EH, a shock is formed between the 
two hot gas halos. The cold disk gas is compressed by 
the hot halo and produces more stars than the other 
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two runs before the closest approach. The dissipation 
of the cold gas occurs much more effectively than the 
accretion of the halo gas onto the disk. When LH swings 
around EH near the closest approach, the inner trailing 
shock propagates toward the center of EH before the tidal 
bridge develops out of the gas disk, which has already 
been largely dissipated, toward EH and blocks its motion 
to the center of EH and ionizes it. Therefore, only the 
old disk stars flow along the tidal bridge, and there is no 
burst of star formation within the cold gas that can be 
accreted to EH. The result of the interaction is mostly 
the old disk star accretion from LH to EH. The SER 
of this run arises somewhat before the closest approach 
because of the shock between the two gas halos. The 
maximum value of the SER after the closest approach 
is much greater than that in the other two runs, but it 
decreases rapidly because of the lack of the cold gas. 

The almost no gas accretion from the LTG to the ETG 
in run EH-LH, unlike in run EH-L, is hinted at in the 
snapshots taken shortly after the closest approach when 
the gas bridge develops. We show in Eigure 11 the gas 
bridges in runs EH-L and EH-LH at t = 5 Gyr more 
closely. In the top panels, the gas particles within the 
disks are displayecQ on a new x—y plane that has an 
X axis along the line connecting the centers of the two 
galaxies and a zero point at the position where the net 
local gravitational acceleration is zero. The gas bridge in 
run EH-L (left panel) develops toward the ETG, well be¬ 
yond the X = 0 line at the time. However, in run EH-LH 
(right panel), the gas bridge does not grow well, reach¬ 
ing behind the x = 0 line. In order to transfer some gas 
to the ETG in run EH-LH, the gas bridge should have 
grown beyond the zero-gravity line before it is blocked 
by the shock. In addition, the pressure of the gas par¬ 
ticles, shown in the bottom-right panel (magenta solid 
line), rises near the end of the bridge. This implies that 
those particles feel the net force (gravity plus pressure 
gradient force) toward the LTG and will not be able to 
flow to the ETG at later times. 

The above results would vary in different situations. 
Eor example, in a case where LH encountered with EH 
more closely (keeping the other orbital parameters un¬ 
changed), the gas tidal bridge might be able to grow 
more strongly to the ETG before being blocked by the 
shock. Then, some gas will manage to flow to the ETG 
despite the shock. The materials transferring through 
the gas and stellar bridges might include star-forming 
gas and young stars in run EH-LH, as in run EH-L. It 
is also possible for the non-star-forming cold gas to flow 
to the ETG for a while through the gas bridge and turn 
into stars there when they become dense enough. The 
gas bridge is expected to vanish more quickly in run EH- 
LH after all, compared with the gas bridge in run EH-L, 
because the shock would still cross the gas bridge and 
ionize it. In the case of EH-L, the SER might start to 
increase before the closest approach, as in run EH-LH, 
when they encountered each other more closely. Einally, 
in run E-L we expect that the disk materials captured by 
the ETG will hardly include star-forming gas and young 

The original disk gas particles that are stripped off or scat¬ 
tered into the halo are not plotted. Those particles that joined the 
halo are easily distinguished on the pgas—T plane, as explained in 
Paper 1. 


stars, even in somewhat closer encounters because there 
is still no compression of the cold gas by shock. The ma¬ 
terials transferring to E through the star and gas bridges 
would also more likely orbit around the center of E and 
rarely settle down at the center, unless all of the model 
and orbital parameters were matched to cancel out the 
transverse component of the velocities of the particles, 
leaving only the radial component. Some cold gas trans¬ 
ferred to E might be able to turn into stars later when it 
becomes dense enough in closer encounters. 

Overall, our simulations show that the cold gas flow 
and the star formation activity depend strongly on 
whether or not the hot gas halo is included in the ETG 
and/or the LTG, whereas the amount of the old disk 
stars accreted to the ETG and their dynamics are al¬ 
most the same in all runs. Therefore, we claim that it 
is necessary to include hot halo gas properly in galaxy 
interaction simulations in order to more accurately study 
the transient phenomena. 
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APPENDIX A 

EVOLUTION OE MODEL LH IN ISOLATION 

In order to check the stability of our galaxy models, 
we have run each of our galaxy models in isolation for 
about 7 Gyr. Eigure 12 presents the distribution of the 
star (top two rows) and gas (bottom two rows) particles 
on the disk at t = 0, I, 4, and 6 Gyr. As time passes, 
some gas particles turn into stars, and spiral patterns 
and a central bar develop in the disk (the bar is seen in 
a more close-up view that is not included in this paper). 
At later times, the gas on the disk becomes dissipated. 
This is primarily because the gas has been continuously 
consumed by star formation, and also some of the disk 
gas has been heated by hot halo gas and joined the halo 
(see the blue points detached from the disk), while gas 
accretion from the halo occurs not very actively. The 
amount of the newly formed stars out of the original halo 
gas particles compared to that out of the original disk 
gas particles is less than I %. Among the star-forming 
gas particles (cyan points), the percentage of the original 
halo gas particles that cooled down onto the disk is less 
than I %, as well. Compared to model L without a gas 
halo, the star formation and the gas dissipation in the 
disk in model LH occur more actively over the course of 
the simulation because of the pressure on the gas disk 
exerted by the gas halo, despite of the gas heating. As 
shown in the figure, the centers of both the star and gas 
disks stay at the origin throughout the simulation. We 
have also checked whether the total energy of the model 
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Figure 12. Evolution of the isolated model LH, which possesses both a gas disk and a gas halo as well as a stellar disk, a bulge, and 
a DM halo. The top two rows display the distribution of the star particles on the disk seen in the x—y plane (first row) and in the x—z 
plane (second row) at t = 0, 1, 4, and 6 Gyr (from the first to the fourth columns, respectively). Orange points represent the original disk 
star particles set from the start of the simulation. Light-orange and green points denote stars formed out of the gas during the simulation 
before t = 4.75 Gyr and since t = 4.75 Gyr, respectively. The bottom two rows show the corresponding views of the gas disk. Blue points 
represent the original disk gas particles with a zero star formation rate. Gyan points are the gas particles having positive values of star 
formation rates. 
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Figure 13. The same snapshots as in Figure 2, but from run E-LH. The pink points in the bottom panels represent the particles originally 
set as the halo gas of the late-type galaxy LH. 
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Figure 14. The local density and the temperature of the gas particles in the high-resolution runs of EH-L (top two rows) and EH-LH 
(bottom two rows) at t = 4 and 5 Gyr. As in Figures 7 and 9, all gas particles within | 2 ;| < 15 kpc are displayed, with the same ranges of 
the color bars. 
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system is conserved while it evolves. The maximum vari¬ 
ation of the total kinetic and potential energies is less 
than 0.9 % and 0.4 % of the initial values, respectively. 
The maximum change of the total internal energy of all 
of the gas particles is less than 7.0 % of the initial value. 

APPENDIX B 

EVOLUTION OE RUN E-LH 

We have considered the encounter between E and LH, 
where only the less massive LTG has a gas halo, for 
comparison and completeness of our study. All of the 
initial orbital and simulation parameters are fixed as in 
the other three runs. In run E-LH, the closest approach 
occurs at t = 4.74 Gyr with a separation between the 
galaxies of 94.5 kpc. Eigure 13 shows the snapshots of 
this run at the same times as those in Eigure 2. The 
distribution of the old disk stars (orange points) is al¬ 
most identical to that in the other three runs. At the 
second snapshot, the newly formed stars since the clos¬ 
est approach (green points) along the bridge appear less 
prominently compared with run EH-L, where the cold 
disk collides directly with the hot halo of EH. In this 
run, many young stars are found in the disk because of 
the compression by the gas halo of LH, as in run EH-LH. 
The SER rises as high as in run EH-LH (see the bottom 
panel in Eigure 10), but it increases abruptly right after 
the closest approach. As shown in the bottom panels 
of Eigure 13, the halo gas of the LTG (pink points) is 
strongly disturbed by the gravitational pull of the more 
massive ETG. The initial disk gas particles that have 
detached from the disk follow the motion of the nearby 
halo gas of the galaxy. 

APPENDIX C 
HIGH-RESOLUTION RUNS 

We have performed our simulations with a higher reso¬ 
lution as well, in order to check particularly whether the 
shock features previously seen in runs EH-L and EH-LH 
were properly resolved. In the high-resolution runs, we 
use the LTG models with a four times greater number 
of particles, keeping all of the other model and orbital 
parameters unchanged. (Eor efficient computing perfor¬ 
mance, we used the same number of particles for the 
ETG models.) Eigure 14 presents the local density and 
the temperature of the gas particles near the orbital plane 


(within ki < 15 kpc) from the high-resolution runs. As 
seen in the figure, the shocks occurring between the gas 
disk and the gas halo in run EH-L and between the two 
gas halos in run EH-LH are almost identical to those in 
Eigures 7 and 9. Because the shocks strongly lead the 
rest of the evolution, the main results of our simulations 
remain consistent in different resolution runs. 
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